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ABSTRACT 

A steady boundary layer flow of Powell- Eyring nanofluid past a stretching sheet with variable thickness in the 
presence of thermal radiation and viscous dissipation is studied numerically. The model is used for the nanofluid 
incorporates the effects of Brownian motion and thermophoresis. The suitable transformations are applied to convert the 
governing partial differential equations into a set of nonlinear coupled ordinary differential equations. Runge-Kutta-based 
shooting technique is employed to yield the numerical solutions for the model. The obtained results for the velocity, 
temperature and concentration are analyzed graphically for several physical parameters. It is found that an increment in 
wall thickness parameter results in decrease of velocity, temperature and concentration profiles. Further, in tabular form 
the numerical values are given for the local skin friction coefficient, local Nusselt number and Sherwood number. A 
remarkable agreement is noticed by comparing the present results with the results reported in the literature as a special 
case. 
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INTRODUCTION 

Viscous boundary layer flow due to a stretching/shrinking sheet is of significant importance due to its vast 
applications. Aerodynamic extrusion of plastic sheets, glass fiber production, paper production, heat treated materials 
travelling between a feed roll and a wind-up roll, cooling of an infinite metallic plate in a cooling bath and 
manufacturing of polymeric sheets are some examples for practical applications of non-Newtonian fluid flow over a 
stretching/shrinking surface. The quality of the final product depends on the rate of heat transfer at the stretching 
surface. This stretching/shrinking may not necessarily be linear. It may be quadratic, power-law, exponential and so 
on. Bilal Asharf et.al. [1] Studied the three dimensional boundary layer flow of Eyrng-Powell nanofluid by 
convectively heated exponentially stretching sheet. Hayat et.al. [2] Investigated the axisymmetric Powell-Eyring 
fluid flow with convective boundary condition. Javed et.at. [3] Analyzed the boundary layer flow of an Eyring - 
Powell non-Newtonian fluid over a stretching sheet. Hayat et.al. [4] Studied the effect of MHD boundary layer flow 
of Powell-Eyring nanofluid over a non-linear stretching sheet with variable thickness. Nazar et.al. [5] Investigated 
the mixed convection boundary layer flow an isothermal horizontal circular cylinder embedded in a porous medium 
filled with a nanofluid for both cases of a heated and cooled cylinder. 

Solar energy is probably the most suitable source of renewable energy that can meet the current energy 
requirements. The energy obtained from nature in the form of solar radiations can be directly transformed into heat 
and electricity. The idea of using small particles to collect solar energy was first investigated by Hunt [6] in the 
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1970s. Researchers concluded that with the addition of nanoparticles in the base fluids, heat transfer and the solar 
collection processes can be improved. Masuda et al. [7] discussed the alteration of thermal conductivity and viscosity by 
dispersing ultra-fine particles in the liquid. Choi and Eastman [8] were the first to introduce the terminology of nanofluids 
when they experimentally discovered an effective way of controlling heat transfer rate using nanoparticles. Buongiorno [9] 
developed the nonhomogeneous equilibrium mathematical model for convective transport of nanofluids. He concluded that 
Brownian motion and thermophoretic diffusion of nanoparticles are the most important mechanisms for the abnormal 
convective heat transfer enhancement. The relevant processes are briefly described in [10-12]. Investigations in the 
nanofluid flows have received remarkable popularity in research community in last couple of decades primarily due to their 
variety of applications in power generation, in transportation where nanofluid may be utilized in vehicles as coolant, shock 
absorber, fuel additives etc., in cooling and heating problems which may involve the use of nanofluids for cooling of 
microchips in computer processors, in improving performance efficiency of refrigerant/air-conditioners etc. and in 
biomedical applications in which magnetic nanoparticles may be used in medicine, cancer therapy and tumor analysis. 
Recently the researchers have proposed the idea of using solar collector based nanofluids for optimal utilization of solar 
energy radiation [13, 14]. Buongiorno and Hu [15] discussed the heat transfer enhancement via nanoparticles for nuclear 
reactor application. Huminic and Huminic [16] showed that use of nanofluids in heat exchangers has advantage in the 
energy efficiency and it leads to better system performance. Ramzan et.al. [17] Discussed effects radiation and MHD on 
Powell-Eyring nanofluid over a stretching cylinder with chemical reaction. Ullah et.al. [18] Studied theoretically the 
influence of thermal radiation and MHD on natural convective flow of casson nanofluid over a non linearly stretching 
sheet. Mustafa et.al. [19] Runge-Kutta fourth -fifth order method with shooting technique to explore the steady boundary 
layer flow of nanofluid past a vertical plate with nonlinear radiation. Numerical solution for steady boundary layer flow of 
dusty fluid over a radiating stretching surface embedded in a thermally stratified porous medium in the presence of uniform 
heat source was studied by Gireesha et.al. [20]. Recently, Ramzan et.al. [21] Analyzed the radiative Jeffery nanofluid with 
convective heat and mass conditions. 

To the best of the author’s knowledge, there seems no existing document about steady boundary layer flow of 
Powell- Eyring nanofluid past a stretching sheet with variable thickness in the presence of thermal radiation and heat 
viscous dissipation. Hence the aim of the present investigation is to examine the steady boundary layer flow of Powell- 
Eyring nanofluid past a stretching sheet with variable thickness in the presence of thermal radiation and viscous 
dissipation. 

MATHEMATICAL FORMULATION 


Radiative powell-Eyring nanofluid over an impermeable nonlinear heated stretching sheet with variable thickness 
is considered. An incompressible fluid is selected electrically conducting. A non-uniform magnetic field 

n —1 

B(x ) = B 0 (x + b ) 2 is imposed transverse to the stretching sheet. Magnetic Reynolds number is chosen small. Induced 

magnetic and electric fields are not accounted. Brownian and thermophoresis in the nanofluid are considered. Newly 
developed condition for mass flux is imposed. The fluid formation is such that the x-axis is presumed along the sheet while 
y-axis is transverse to it. The temperature of the sheet is different from that of the ambient medium. The fluid viscosity is 
assumed to vary with the temperature while the other fluid properties are assumed constant. 


The stretching surface has the nonlinear velocity U w 


U 0 (x + b)" where U 0 


is the reference velocity and b is 
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the dimensional constant. Further it is assumed that sheet at y = A (x + b ) 2 is not flat (where n is the velocity power 
index and A is assumed very small so that the sheet retain adequately thin). We also noticed that for n = l the problem 
reduces to a flat sheet. The governing expressions for considered flow are given by: 

du dv 


- + — = 0 

c)x dy 


du du 

U lh +X %- 


dT dT 
u -1-v- 
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With boundary conditions 

u=U w (x) = U 0 (x + b) n , v = 0, T = T w C = C W at y = A(x + b) 
u —y 0, v —y 0, T —> T^ C —\ C x as y —y °° 


1—77 

T 


(5) 


(6) 


Here u and v are the corresponding velocity components parallel to x- and y-directions respectively, // the 

P 0 

dynamic viscosity, V = — designates the kinematic viscosity, p the fluid density, d and p are the material liquid 

p 


parameters of Powell-Eyring model, k the thermal conductivity, a 


f k A 


nf 


M 


the thermal diffusivity of liquid, T the 


f J 


ratio of the heat capacity of liquid of the nanoparticles material to the effective heat capacity of the base fluid, D B 
indicates the Brownian diffusion coefficient, D T represents the thermophoretic diffusion, T the temperature of the fluid, C 
the nanoparticles concentration, T and T x are the sheet and ambient fluid temperatures and C, t the ambient fluid 
nanoparticles concentration. 
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Figure 1: Physical Model and Coordinate System 

Using the Rossel and approximations, the radiative heat flux is given by 

4(7* dT 4 

q r = --— 

3K* dz 


(7) 


Where, (7 * and K * are, respectively, the Stephan - Boltzman constant and the mean absorption coefficient. 
Assume that the temperature differences within the flow are small such that T 4 in a Taylor series about T x and neglecting 
higher order terms. We get 

T 4 = 4TT - 3T 4 (8) 


Substituting eq.’s (7) and (8) in eq. (3). We have 
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Transformations are expressed as follows: 
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Incompressibility condition is satisfied identically and Eqs.(2), (4)-(6) and (9) take the following forms 
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<E> A ’ + Pr LeF<t>' + 
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In the above expressions primes designate differentiation with respect to 7J . Where a = A 


(13) 
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wall thickness parameter and Tj = a = A 


'n + O 


v z y 


u o 

- shows the plate surface. Upon using 
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( 20 ) 


where N and A are the fluid parameters, M represents the magnetic parameter, Nb the Brownian motion 
parameter , Pr the Prfzandtl number, Nt the thermophoresis parameter, Le presents Lewis number, Nr the radiation 
parameter, Ec the Eckert number and prime indicates differentiation with respect to f. The non-dimensional parameters 
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Surface drag coefficient and surface heat transfer are expressed as follows: 
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C = Nu - (x+b)q ' 
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In dimensionless form we have [37, 38]: 

C, Re> = J2(n + l)f(1 + W)/'(0)-<^3)M(/-(o)) 3 ' 

V 1 $ ) 


(24) 


AfeRe; K =-,p4V(0) 


(25) 


Where Re designates the local Reynolds number presented by 


U w {x + b) 


SOLUTION OF THE PROBLEM 


The governing equations for the present problem are transformed into a set of coupled nonlinear differential 
equations by applying similarity transformation. The equations (16) - (18) together with the boundary conditions (19) - 
(20) are integrated numerically by using Runge-Kutta-Gill method along with the shooting technique. This method is 
concisely outlined as below: 
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The boundary conditions are transformed as follows: 


>’, (0) = a 


1-/7 
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, > 2 (0) = 1, > 4 (0) = 1, > 6 (0) = 1 


\i + nj 
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(27) 


In order to carry out the step by step integration for the equations (16) - (18), Gill’s procedures have been used 
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(Ralston and Wilf [23]). To start the integration it is necessary to provide all the values of y x , y 2 , y 3 , y 4 , y 5 , y 6 at // = 0 
from which point, the forward integration has been carried out but from the boundary conditions it is seen that the values of 
y 3 , y 4 , y n are not known. So, we are to provide such values of y 3 , y 4 , y 7 along with the known values of the other function 
at Tj — 0 as would satisfy the boundary conditions as tj -4 °° to a prescribed accuracy after step by step integrations are 
performed. Since the values of y 3 , y 4 , y 7 which are supplied are merely rough values, some corrections have to be made in 
these values in order that the boundary conditions to 77 —> are satisfied. These corrections in the values of y 3 , y 4 , y 7 are 
taken care of by a self-iterative procedure which can for convenience be called corrective procedure. 

RESULTS AND DISCUSSIONS 

In this section, we have presented the physical behavior of the velocity, temperature and concentration profiles for 
different physical parameters such as magnetic parameter M, fluid parameters A, N , wall thickness parameter CC, Prandtl 
number Pr, Eckert number Ec, thermophoresis parameter Nt, Brownian motion parameter Nb and thermal radiation 
parameter R. we are solving the system of coupled nonlinear equations by using well known method Runge-Kutta-Gill 
method along with the shooting technique. 

Figure 2 portrays velocity field for varied values of magnetic parameter M and fluid parameters N. It is observed 
that velocity profiles reduces with an increase in magnetic parameter M, physically, presented results occurs in the 
conductivity fluid and create a resistive type force i.e. Lorentz force on the fluid in the boundary layer , which fall 011 the 
velocity of the fluid . In the presence of fluid parameter, it is noteworthy that the velocity profile depicts an increment with 
increasing N. The behavior of M and N on the temperature and concentration is shown in figs.3 and 4 respectively. 
Temperature and nano concentration profiles show increasing behaviour for higher values of M, whereas decreasing 
behaviour for N. An increment in M produces higher Lorentz force (resistive force) which has the characteristic to convert 
some thermal energy into heat energy. Therefore the temperature and concentration profile increases. Furthermore, the 
opposite behaviour was observed for the effect of N . Physically, it justifies that higher value of N produce a reduction in 
the fluid viscosity. Therefore, less heat is produced due to frictional force. Hence, the temperature and concentration 
distribution decreases. 

Figure 5 displays that the wall thickness parameter Ot decreases velocity profile of the fluid physically, increasing 
the value of CC decrease the flow velocity, because under the variable wall thickness, not all the pulling force of the 
stretching sheet can be transmitted to the fluid causing a decrease for the friction between the fluid layers. Effects of the 
thermophoresis parameter Nt and Brownian motion parameter Nb on temperature and nano concentration distributions 
are shown in figures 6 and 7. It is comprehended that temperature and concentration fields are increasing functions of 
thermophoretic parameter. In fact, small particles are pulled towards cold surface from hot one. Eventually, temperature of 
the fluid increases which results in higher temperature and nano concentration profiles. From these figures, it also analyzed 
that temperature and concentration distributions show increasing and decreasing behaviour for higher values of Brownian 
motion parameter with upsurge in values of Brownian motion parameter an increase in fluid particles random motion is 
observed that eventually produce more heat. So, we observe increasing temperature field and decreasing concentration 
distribution. Figure 8 shows that the effects of thermal radiation parameter and Eckert number on temperature profiles. R = 
0 shows that no radiation on the fluid layer and Ec = 0 represents the no viscous dissipation on the fluid layer. The 
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temperature profile increases with an increase in the thermal radiation and Eckert number. 

Furthermore, from the figure 9 and 10, it is observed that an increase in the Prandtl number results in decreasing 
the heat and mass transfer profiles. The reason is that increasing values of Prandtl number is equivalent to decreasing the 
thermal conductivities, and therefore heat is able to diffuse away from the heated sheet more rapidly. Hence, in the case of 
increasing Prandtl number the boundary layer is thinner and the heat transfer is reduced. The effect of Lewis number on 
concentration profile is presented in figure 11. It is observed that concentration profiles decreases with increase in Le. 
Figures 12 and 13 depicted for the effects of thermophoresis parameter and Brownian motion parameter on local Nusselt 
number and local Sherwood number, respectively. From these figures, it is noticed that local Nusselt number decreases 
with an increase in Nt and Nb. The local Sherwood number decreases with an increase in Nt, whereas it increases with an 
increase in Nb. 

To verify the accuracy of the numerical results, we compared our results with those reported by Akbar et.al. [24] 
and Malik et.al. [25] (for regular fluids) as shown in table 1. The results are very good agreement, thus lending confidence 
to the accuracy of the present results. 

CONCLUSIONS 

We inspected the feature of radiative powell-Eyring nano fluid past a stretching sheet with variable thickness and 
viscous dissipation. An increment in wall thickness parameter results in decrease of velocity, temperature and 
concentration profiles. Enhancement in the fluid temperature is observed for higher thermophoresis parameter. 
Temperature profile increases for increasing values of radiation and viscous dissipation. Magnitude of both local Nusselt 
number and Sherwood numbers are decreasing functions of Nt. 



Figure 2: Velocity Profiles for Various Values of M & N. 
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Figure 3: Temperature Profile for Various Values of M & N. 



T1 

Figure 4: Concentration Profiies for Various Values of M & N. 



T| 

Figure 5: Velocity Profile for Various Values of (X 
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T1 


Figure 6: Temperature Profiles for Various Values of Nt & Nb. 



Figure 7: Concentration Profiles for Various Values of Nt & Nb. 



n 

Figure 8: Temperature Profile for Various Values of Ec & R. 
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Figure 9: Temperature Profiles for Various Values of Pr. 



Figure 10: Concentration Profiles for Various Values of Pr. 



Figure 11: Concentration Profiles for Various Values of Le. 
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Figure 12: Local Nusselt number for various values of Nt & Nb. 



Figure 13: Local Sherwood Number for Various Values of Nt & Nb. 


Table 1: Comparison of Skin-Friction Coefficient C ; Re' 7 with the Available 
Results in Literature for different Values of M when n = X — Ot — 0 



C f Re 

1/2 


M 

Present Study 

Akbar et al. [24] 

Malik et al.[25] 

0 

1 

1 

1 

0.5 

-1.118034 

-1.11803 

-1.11802 

1 

-1.414214 

-1.41421 

-1.41419 

5 

-2.449490 

-2.44949 

-2.44945 

10 

-3.316625 

-3.31663 

-3.31657 

100 

-10.049874 

-10.04988 

-10.04981 

500 

-22.383029 

-22.38303 

-22.38294 

1000 

-31.638584 

-31.63859 

-31.63851 
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